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Abstract. We report optical extinction properties of dust for a sample of 26 early-type galaxies based on the analysis of 
their multicolour CCD observations. The wavelength dependence of dust extinction for these galaxies is determined and the 
extinction curves are found to run parallel to the Galactic extinction curve, which implies that the properties of dust in the 
extragalactic environment are quite similar to those of the Milky Way. For the sample galaxies, value of the parameter Ry, 
the ratio of total extinction in V band to selective extinction in B & V bands, lies in the range 2.03 - 3.46 with an average of 
3.02, compared to its canonical value of 3.1 for the Milky Way. A dependence of Ry on dust morphology of the host galaxy 
is also noticed in the sense that galaxies with a well defined dust lane show tendency to have smaller Ry values compared 
to the galaxies with disturbed dust morphology. The dust content of these galaxies estimated using total optical extinction is 
found to lie in the range 10* to 10* Mq, an order of magnitude smaller than those derived from IRAS flux densities, indicating 
that a significant fraction of dust intermixed with stars remains undetected by the optical method. We examine the relationship 
between dust mass derived from IRAS flux and the X-ray luminosity of the host galaxies.The issue of the origin of dust in 
early-type galaxies is also discussed. 
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1. Introduction 

Recent imaging surveys using both ground and space-based 
telescopes across the electromagnetic spectrum have amply 
demonstrated that early-type galaxies, ellipticals (E) and lentic- 
ulars (SO) contain complex, multi-phase interstellar medium 
(ISM). In particular, a large fraction (~ 50 - 80%) of E-SO 
galaxies, at least in our local universe, are now known to pos- 
sess dust features in a variety of morphological forms, as re- 
vealed from high quality optical imaging, thanks to the avail- 
ability of sensitive array detectors (Goudfrooij et al. 1994b,c; 
van Dokkum & Franx 1995; Tran et al. 2001; and refer- 
ences therein). Dust attenuates the radiation at optical wave- 
lengths and re-radiates at infra-red (IR) wavelengths, espe- 
cially in far-infrared (FIR) bands. FIR emission from early-type 
galaxies detected using the Infra Red Astronomical Satellite 
(IRAS) (Knapp et al. 1989), and now with the Infrared Space 
Observatory (ISO) and Spitzer Space Telescope is attributed to 
thermal emission from dust having a wide variety of tempera- 
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tures (Leeuw et al. 2004; Temi et al. 2004; Xilouris et al. 2004). 
The dust masses derived from ISO data are at least an order 
of magnitude larger than those estimated using IRAS fluxes, 
which in turn have been found to be on average an order of 
magnitude higher than the dust mass estimated using optical 
extinction. 

The study of physical properties of dust such as extinction, 
reddening and polarization help in deriving important informa- 
tion about the total dust content, size of the dust grains respon- 
sible for extinction and its variation with environment, metal- 
licity, star formation history and redshift of the host galaxy. 
Dust properties in galaxies at high redshift are noticeably dif- 
ferent from those found at later cosmic times, because the dust 
grains are smaller in size due to their different formation his- 
tory (through type II Supernovae) and also due to the short 
time available for accreting heavy atoms and coagulation with 
other grains (Maiolino et al. 2004). Thus, detailed study of the 
dust properties in the extragalactic environment provides useful 
clues not only for understanding the origin and fate of dust in 
external galaxies, but also for the subsequent evolution of host 
galaxies. 
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Our knowledge of dust properties largely relies on the inter- 
action of dust particles with the electromagnetic radiation i.e., 
on attenuation and scattering of starlight, collectively known as 
extinction, and re-radiation by dust at longer wavelengths. The 
spectral dependence of extinction, termed as extinction curve, 
is found to depend strongly on the composition, structure and 
size distribution of the dust grains. Therefore, the first step in 
understanding the physical properties of dust in the extragalac- 
tic environment is to derive the extinction curves for the sample 
galaxies. For our own Galaxy, the /l-dependence of extinction 
is derived by comparing the spectral energy distribution (SED) 
of a pair of stars of identical spectral and luminosity class with 
and without dust in front of them (Massa et al. 1983). Any dif- 
ference in the measured magnitude and colour of these stars 
is attributed to the dust extinction. The standard extinction law 
thus derived is found to be uniformly applicable in our Galaxy 
from optical to near-infrared wavelengths. Over the UV to far- 
IR wavelength range the standard extinction law is character- 
ized by a single parameter Ry, the ratio of total extinction Ay in 
V band to the selective extinction in Z? & V bands i.e., E{B - V), 
with Rv = 3.1 (Savage & Mathis 1979; Rieke & Lebofsky 
1985; Mathis 1990). However, within the Milky Way, the value 
of Ry is found to vary between 2. 1 to 5.6 depending on the line 
of sight (Valencic et al. 2004). Similarly, the extinction curves 
have also been derived for a few neighbouring spiral galaxies 
using this method and the Ry value is found to vary consider- 
ably for these galaxies (Brosch 1988), and even within a galaxy 
from one location to another 

In the case of external galaxies, as individual stars can- 
not be resolved, the above method is not applicable; instead, 
a variety of indirect methods have been proposed to estab- 
lish the /l-dependence of the dust extinction, and to estimate 
Ry values. In one of the widely used methods, comparison 
of light distribution of the original galaxy with its dust-free 
smooth model gives an estimate of the extinction caused by 
dust present in the original galaxy. The extinction law for ex- 
ternal galaxies can be deduced by performing this exercise at 
different wavelengths. The absence of spiral arms, H ii regions 
and other inhomogeneities results in a fairly smooth light dis- 
tribution in early-type galaxies. This allows one to easily and 
accurately construct the dust-free model of the galaxy required 
to study dust extinction properties. Despite this rather conve- 
nient situation, dust extinction properties have been investi- 
gated only for a handful of galaxies, most notably by Brosch 
& Loinger (1991), Brosch & Almoznino (1997), Goudfrooij et 
al. (1994b,c) and extended to a few more individual galaxies by 
Sahu et al. (1998), Dewangan et al. (1999), Falco et al. (1999), 
Keel & White (2001), Motta et al. (2002) etc. In all these stud- 
ies it has been observed that the optical extinction curves in the 
extragalactic environment closely resemble that of the Milky 
Way, with Ry values comparable to the canonical value of 3.1. 
Dust being an important component of the ISM, it is essential 
to extend this kind of analysis to a larger sample of dusty early- 
type galaxies. 

It has been shown that the morphology of dust closely 
matches that of the ionized gas in a large fraction of galaxies (~ 
50 - 80 %) (Goudfrooij et al. 1994b, Ferrari et al. 1999), and in 
some cases with the X-ray emitting region too (Goudfrooij & 



Trinchieri 1998), pointing to a possible physical connection be- 
tween hot, warm and cold phases of ISM in early-type galaxies. 
In the multiphase ISM, dust grains can act as an efficient agent 
in transporting heat from hot gas to cold gas giving rise to the 
observed warm phase of ISM (Trinchieri et al. 1997). Again a 
detailed investigation of dust properties is called for 

We have an ongoing program of detailed surface photomet- 
ric study of a large sample of early-type galaxies containing 
dust to investigate dust properties in the extragalactic environ- 
ment and compare them with those of the Milky Way. This 
paper reports on dust properties in a sample of 26 early-type 
galaxies, based on their deep, broad band optical {BVRI) imag- 
ing observations. §2 describes the sample selection, observa- 
tions and preliminary reduction of the acquired data, §3 de- 
scribes the properties of dust for the sample galaxies, in §4 re- 
sults obtained are discussed. Our results are summarized in §5. 
Through out this paper we assume Ho - 50 km s ' Mpc 

2. Observations and data reduction 

We have carried out multiband optical imaging observations 
of 26 dusty early-type galaxies (HE and 15 SO) as a part 
of our ongoing program studying dust properties in a large 
sample of early-type galaxies taken from Ebneter & Balick 
(1985), Veron-Cetty & Veron (1988), Knapp et al. (1989), van 
Dokkum & Franx (1995). The objects were chosen depending 
on the availability of observing nights and weather conditions, 
and as such no strict criterion has been applied for the sam- 
ple selection. Table[T]gives the global parameters such as coor- 
dinates, morphological type, heliocentric velocity, luminosity, 
size and environment of the target galaxies. 

Deep CCD images of the program galaxies were obtained 
using various observing facilities available in India, during 
December 1998 to August 2003. Table |2] gives details of the 
instruments used during different observing runs and Table [3] 
gives the log of observations. Except for the lAO observing 
run, where Bessel U,B,V,R,I filters were used, observations 
were made in Johnson B, V and Cousins R, I filters. Generally, 
exposure times were adjusted so as to achieve roughly equal 
signal-to-noise (S/N) ratio for a galaxy in different bands. Apart 
from the object frames, several calibration frames such as bias, 
twilight sky flats etc. were also taken in each observing run. 
For photometric calibration of the data, open cluster M67 and 
standard stars from Landolt's list (Landolt 1992) were observed 
during photometric nights. 

Standard preprocessing steps such as bias subtraction and 
flat fielding were done using the standard tasks available within 
IRAlfl Multiple frames taken in each filter were geometrically 
aligned to an accuracy better than one tenth of a pixel by mea- 
suring centroids of several common stars in the galaxy frames, 
and were then combined to improve the S/N ratio. This also en- 
abled easy removal of cosmic ray events. The cosmic ray hits 
that were left after combining the frames were further elimi- 
nated using the cosmicrays task in IRAF. Sky background in 

' IRAF is distributed by the National Optical Astronomy 
Observatories (NOAO), which is operated by the Association of 
Universities, Inc. (AURA) under co-operative agreement with the 
National Science Foundation. 
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Object 


RA 


DEC 


Morph. 


oil 


Vhelio 


size 


CfA Group 




(JzOOO.O) 


(JzUOO.O) 


KC3(KiA) 




(km/s) 


(arcmin) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 524 


01:24:47 


09:32:21 


SO(SO/a) 


11.17 


2509 


2.8x2.8 


GH13 (8), dominant in group 


NGC 984 


02:34:43 


23:24:46 


SO 


13.35 


4435 


3.0x2.0 


- 


NGC 1172 


03:01:36 


-14:50:13 


EI(SOi) 


12.56 


1502 


2.3x1.8 


- 


NGC 1439 


03:44:50 


-21:55:21 


El(£l) 


12.18 


1568 


2.5x2.3 


HG32 (21), Eridanus group 


NGC 2128 


06:04:34 


57:37:40 


SO 


12.66 


3226 


1.5x1.1 


- 


NGC 2534 


08:12:54 


55:40:24 


El 


13.38 


3732 


1.4x1.2 


Isolated galaxy 


NGC 2563 


08:20:35 


21:04:09 


SO 


13.01 


4592 


2.1x1.5 


In dense part of Cancer cluster 


NGC 2672 


08:49:22 


19:04:30 


EI (£2) 


12.12 


4164 


3.0x2.8 


Brightest in group, pair with N2673 


NGC 2693 


08:56:59 


51:20:51 


E3(£2) 


12.61 


4923 


2.6x1.8 


GH39 (4) 


NGC 2907 


09:31:36 


-16:44:09 


SO/a(SO) 


12.35 


1892 


1.8x1.1 


Brightest in group, several dEs companion 


NGC 29 11 


09:33:46 


10:09:09 


SO(SO) 


12.25 


3062 


4.1x3.2 


GH47 (6) 


NGC 3078 


09:58:24 


-26:55:34 


E2(£3) 


11.94 


2283 


2.5x2.1 


HG29(4), in cluster 


NGC 3489 


11:00:18 


13:54:08 


SO(SO/a) 


11.15 


613 


3.5x2.0 


GH68 (23)/HG56 (30), Leo grp 




1 1 '07' 18 


-19'28'21 


SO 


13 03 


3672 


2 6x1 4 


In cluster, interacting pair 


NGC 3585 


11:13:17 


-26:45:20 


E6(£/50) 


10.64 


1206 


4.7x2.6 




NGC 3599 


11:15:27 


18:06:46 


SO 


12.70 


713 


2.7x2.1 


GH77 (13) 


NGC 3665 


11:24:43 


38:45:47 


S0(5'0) 


11.69 


2072 


2.5x2.0 


GH79 (4), pair with N3658 


NGC 3923 


11:51:02 


-28:48:23 


E4(£/50) 


10.62 


1487 


5.9x3.9 


HG28 (4), at least 7 dEs 


NGC 4459 


12:29:00 


13:58:46 


S0(5'0) 


11.21 


1154 


3.5x2.7 


GH106 (248) 


NGC 5363 


13:56:07 


05:15:19 


10(50) 




1113 


4.1x2.6 


HG55 (3), pair with N5364 


NGC 5485 


14:07:11 


55:00:08 


SO(SO) 


12.31 


2107 


2.3x1.9 


HG77 (4), pair with N5484, N5486 


NGC 5525 


14:15:39 


14:17:02 


SO 


13.78 


5569 


1.4x0.9 


GH134 (8) 


NGC 5898 


15:18:13 


-24:05:49 


EO(SO) 


11.92 


2153 


2.2x2.0 


In pair with NGC 5903, 
several dEs, dSOs companion 


NGC 5903 


15:18:36 


-24:05:06 


E2(£/50) 


12.20 


2565 


2.7x2.1 


In pair with NGC 5898, 
several dEs, dSOs companion 


NGC 7432 


22:58:06 


13:08:00 


E 


14.41 


7615 


1.5x1.2 


In pair 


NGC 7722 


23:38:41 


15:57:15 


SO 


13.19 


4197 


1.7x1.4 


In pair with a spiral 



Notes: Column (2) and (3) list galaxy co-ordinates, Column (4) lists moiphological classification of the galaxies, RSA type is given in the parenthesis. 
The total, coiTected blue luminosity of the program galaxies are listed in Column (5), while Column (6) lists heliocentric velocity and Column (7) 
optical size of the galaxies, all taken from RC3 (de Vaucoulers et al. 1991). Column (8) lists CfA group membership HG or GH listed by Huchra & 
Geller (1982; HG82) and Geller & Huchra (1983; GH83), number in the parenthesis represents number of galaxies in a given group. 



Table 2. Details of telescopes and instrumentation 



Observing Run 


Dec.98 


Dec.99 


Mar.2k 


Mar.Ol 


Apr.Ol 


May 03 


Aug.03 


Observatory 


ARIES 


ARIES 


VBO 


VBO 


VBO 


lAO 


lAO 


Telescope size(in m) 


1.04 


1.04 


2.34 


2.34 


2.34 


2.0 


2.0 


CCD Type 


TEK 


TEK 


Phot 


Phot 


Phot 


Phot 


Phot 


Format (# pixels) 


Ux \k 


2itx2*; 


litx U 


Ik X \k 


\kxlk 


Ux 1* 


lkx\k 


Binning 


2x2 


2x2 


1 X 1 


1 X 1 


1 X 1 


I x I 


1 X 1 


Scale ("A)in) 


0."73 


0."73 


0."72 


0."72 


0."72 


0."285 


0."285 


Field of view 


13' X 13' 


13' X 13' 


10' X 10' 


10' X 10' 


10' X 10' 


5' x5' 


5' x5' 



NoiEs: Telescope - 1.04 = 1.04m Sampurnanand Telescope of Aryabhatta Research Institute of Observational Sciences 
(ARIES), Naini Tal; 2.34 = 2.34m Vainu Bappu Telescope of V.B. Observatory (VBO), Kavalur; 2.0 ± 2.0m Himalayan 
Chandra Telescope of Indian Astronomical Observatory (lAO), Mt. Saraswati, Hanle. 
CCD type - TEK: Tektronics, Phot : Photometries. 



the galaxy frame was estimated using the box method (see, 
e.g., Sahu et al. 1998); median of a 5x5 pixels box at vari- 
ous locations in the frame, generally away from the galaxy and 
not affected by the stars, was estimated and its mean value was 
taken as a measure of the sky background which was then sub- 
tracted from the respective galaxy frame in the corresponding 
band. As the total field coverage of the CCDs used was large 
compared to the optical size of the program galaxies, the box 
method for sky estimation was found suitable for the present 
study. Cleaned, sky subtracted B,V,R,I images of individual 
galaxies were convolved with a Gaussian function to match the 
seeing of the best frame with that of the worst frame to con- 
struct colour-index images of the galaxies. 

3. Properties of dust: the method and results 

Even though the program galaxies were known to contain 
dust features, we re-examined them using a variety of image 



processing techniques like quotient image, unsharp masking, 
colour-index image, etc., to confirm the presence of dust fea- 
tures and their spatial distribution. (B - V) and (B - R) colour 
index images were constructed using geometrically aligned, 
seeing matched direct images. {B - V) colour index images 
of some of the dust lane galaxies are shown in Fig. \T\ where 
brighter shades represent dust occupied redder regions. Our 
analysis confirms the presence of dust in all the 26 galaxies 
studied here. Dust is present in variety of forms; five galaxies 
have multiple dust lanes parallel to the major axis, twelve ob- 
jects show a well-defined dust lane aligned either along major 
or minor axes, three galaxies show dust rings or arcs, while oth- 
ers show nuclear dust patches. NGC 2907 and NGC 7722 have 
at least four extended dust lanes running parallel to the optical 
major axis. NGC 5363 has two lanes, the inner one is short and 
aligned along minor axis, while the other one is extended and 
lies parallel to the major axis of the galaxy. 
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Table 3. Observing log 
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Galaxy 




Exposure time (sec) 






Seeing 


(") 




Notes 




B 


V 


R 


I 


B 


V 


R 


1 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


NGC 524 


900(4) 


900(3) 


400(3) 


400(3) 


2.65 


2.41 


2.14 


2.11 


#1 


NGC 984 


900(4) 


600(4) 


400(3)+200(l) 


300(4) 


3.01 


2.72 


2.71 


2.91 


#2 


NGC 1172 


900(4) 


900(3) 


400(3) 


400(3) 


2.66 


2.55 


2.46 


2.61 


#2 


NGC 1439 


600(4) 


600(3) 


200(4) 


200(3) 


2.55 


1.93 


1.86 


1.95 


#2 


NGC 2128 


900(4) 


900(3) 


400(2) 


400(2) 


2.81 


2.46 


2.39 


2.09 


#1 


NGC 2534 


900(4) 


900(3) 


600(3) 


600(3) 


2.51 


2.25 


1.79 


2.06 


#1 


NGC 2563 


- 


600(4) 


400(3) 


400(3) 


- 


4.98 


4.52 


4.32 


#3 


NGC 2672 


900(4) 


600(3) 


300(3)+200(1) 


300(3) 


1.90 


1.87 


2.08 


1.88 


#2 


NGC 2693 


900(4) 


600(3) 


360(4) 


300(4) 


3.30 


3.33 


3.07 


2.96 


#4 


NGC 2907 


900(4) 


600(3) 


300(4) 


300(3) 


3.46 


2.61 


2.57 


2.46 


#2 


NGC 2911 


1200(3)+600(1) 


600(4) 


300(5) 


300(4) 


2.62 


2.30 


2.00 


1.86 


#1 








300(2) 




2 87 




3 29 




#4 


NGC 3489 


500(3) 


400(3) 


400(1) 


360(1) 


2.97 


2.75 


2.77 


2.87 


#2 


NGC 3497 




700(3) 




420(2) 




2.38 




2.41 


#7 


NGC 3585 


600(4) 


600(3) 


300(2) 




3.16 


3.07 


2.90 




#3 


NGC 3599 


600(3) 


300(3) 


300(3) 


200(3) 


2.59 


2.54 


2.47 


2.21 


#2 


NGC 3665 


600(3) 


400(3) 


300(2) 




2.11 


1.96 


1.82 




#6 


NGC 3923 


900(4) 


600(4) 


300(4) 




3.24 


3.60 


3.19 




#3 


NGC 4459 


900(4) 


600(3) 


360(4) 




2.50 


2.41 


2.20 




#5 


NGC 5363 


600(4) 


600(3) 


300(4) 


300(2) 


2.12 


1.96 


2.42 


1.82 


#6 


NGC 5485 


600(6) 


600(3) 


300(4) 


300(6) 


2.11 


1.96 


1.54 


1.52 


#7 


NGC 5525 


600(6) 


500(4)+280(l) 






2.01 


2.11 






#7 


NGC 5898 


800(3) 


400(3) 


300(2) 


300(2) 


3.25 


2.61 


3.06 


2.89 


#3 


NGC 5903 


800(3) 


400(3) 


300(2) 


300(2) 


3.25 


2.61 


3.06 


2.89 


#3 


NGC 7432 


600(6) 


600(3) 


300(6) 


300(6) 


1.41 


1.37 


1.25 


1.48 


#6 


NGC 7722 


600(5) 


600(3) 


400(3) 


400(3) 


1.45 


1.75 


1.30 


1.61 


#4 



Notes to Table [3] Column (1) lists sample galaxies. Columns {2)-{5) list exposure time in seconds for B, V, R, & I bands, respectively, 
number in parenthesis represents the number of exposures acquired in a given band. Seeing (in arcsec) in different bands are listed in columns 
(6)-(9). Column (10) lists observing run and instruments used; #1 - ARIES, 1998; #2 - ARIES, 1999; #3 = VBO, March 2000; #4 = VBO, 
March 2001; #5 = VBO, April 2001; #6 = lAO, May 2003; #7 - lAO, August 2003 (see Table [|). 



NGCSDOT 




NGC 348D 


NGC 34B7 




NGC .1665 


^^^^^^^^^^^^^^^^^^^ 
















NGC ,i3K! 




MiC 3185 


NGC i)o23 




NGC 7722 



Fig. 1. (B - V) colour index maps of some of the prominent dust lane galaxies; brighter shade represents dust occupied regions 



3.1. Extinction maps 

To estimate the effect of dust extinction, we compare the light 
distribution in the original galaxy image with its dust free 
model. As early-type galaxies have a fairly smooth and sym- 
metric stellar light distribution with respect to the nucleus, 
one can easily construct its smooth, dust free model by fit- 
ting ellipses to the isophotes of observed image. This method 
has been used by a number of researchers (Brosch & Loinger 
1991; Goudfrooij et al. 1994c; Sahu et al. 1998, and references 
therein) to study extinction properties of dust in the extragalac- 
tic environment. 

We have fitted ellipses to the isophotes of the observed 
galaxy images using the ellipse fitting routine available in the 



STSDAS0 package, which is based on a procedure described 
by Jedrzejewski (1987). 

Starting with trial values of eUipticity, position angle, and 
centre coordinates of the galaxy image, an ellipse was fitted to 
the isophote at a given semi-major axis length after masking all 
the obvious regions occupied by foreground stars and interact- 
ing galaxies, which were ignored during the ellipse fit. The fit- 
ting was continued by incrementing the semi-major axis length 
by 10% until signal reaches 3o" level of the background, and 
continued inward until the centre of the galaxy. A model im- 
age constructed using the best fit ellipses was then subtracted 

- STSDAS is distributed by the Space Telescope Science Institute, 
operated by AURA, Inc., under NASA contract NAS 5-26555. 
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from the original galaxy image and its residual image was gen- 
erated. The residual image was examined and regions occupied 
by dust and other hidden features were flagged and rejected 
in the next run of ellipse fitting. As the majority of galaxies 
from this sample contain complex dust lanes or patches pass- 
ing through their centres, accurate estimation of centre coordi- 
nates was not possible. To minimize the error in this estimate, 
isophote fitting was first carried out for R or I band images, 
which are least affected by dust extinction among the avail- 
able bands and centre coordinates were determined by aver- 
aging those of the best fitted ellipses. Centre coordinates thus 
determined were kept fixed in the second run and the above de- 
scribed ellipse-fitting procedure was repeated. The same centre 



coordinates were used to fit ellipses to the isophotes in other 
bands. 

The "dust free" model of the program galaxies thus gen- 
erated were used to construct "extinction maps" in magnitude 
scale using the relation 



-2.5 X log 



h 



i, model 



(1) 



where represents the amount of extinction in a particular 
band A (= B,V,R,I) while Ii,obs and I^,modei represent ADU 
count levels in the original and model galaxies, respectively. 

Extinction maps for the sample galaxies are shown in 
Fig. |2] where brighter features represent regions of higher op- 
tical depth associated with the dust extinction. Asymmetries 
seen in the dust features of NGC 5903 and NGC 7432 are due 
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to the masking of foreground stars present near the centre of 
galaxies. 

3.2. Extinction curves 

The next step involves the quantification of total extinction 
in each band and deriving the extinction curve. For this pur- 
pose, masks were set on the regions occupied by dust in extinc- 
tion maps generated above, and numerical values of (with 
A-B,V,R&r) were extracted as the mean extinction within a 
square box of 5 x 5 pixels (comparable to the size of the seeing 
disk) and the box was slided over the dust-occupied region in 
each galaxy. To avoid seeing effects, we excluded nuclear re- 
gions (radius < 5") of the program galaxies from analysis. The 
values of extinctions A^(x,y) thus measured at different loca- 
tions in individual bands were used to estimate the local values 
of selective extinction or colour excess E{A - V) = - Ay as 
a function of position across the dust-occupied region. 

A linear regression fit was performed between various local 
values of total extinction {AB,Av,AR,Af) and the slopes of the 
best fits were assigned to be the average slope of A ^ versus Ay 
(where x,y - B,V,R&I;x + y) and the reciprocal slope of A,, 
versus A V (Goudfrooij et al. 1994b; Sahuet al. 1998; Dewangan 
et al. 1999). Likewise, slopes of the fitted lines of the regression 
for A,i {A = B,V,R& I) and E(B - V) were also derived. The 
best fitting slopes were used to derive R,\ ^= -g(§iYj\ for the 
dust occupied regions in the program galaxies and are listed 
in Table |4] along with their associated errors. R^ values for the 
Milky Way taken from Rieke & Lebofsky (1985) are also listed 
in the table for comparison. 

The extinction curves for the program galaxies are given 
in Fig. [3] along with that for the Milky Way. These figures 
demonstrate that the extinction curves for the majority of galax- 
ies run parallel to that of the Galactic curve, except for a few 
cases (NGC 1439, 2534, 2693, 3585, 5363) for which we see 
"concave" extinction curves. The fact that, on average, extinc- 
tion curves for the sample galaxies are parallel to that of our 
Galaxy implies that the dust extinction properties in the extra- 
galactic environment are similar to those of the Milky Way. 
However, Ry for program galaxies is found to be different from 
the canonical value of 3.1 for the Milky Way, as evident from 
Tables 

Assuming that the chemical composition of dust grains in 
the extragalactic environment is similar to that of our Galaxy, 
smaller (larger) Ry values imply that the dust grains responsi- 
ble for interstellar reddening are smaller (larger) than the grains 
in our Galaxy. This allows one to compute the relative grain 
size in the sample galaxies by making use of available mod- 
els of chemical composition and shape of the dust grains. Out 
of several "dust models" a two-component model consisting of 
spherical graphite and silicate grains, with an adequate mix- 
ture of sizes (Mathis et al. 1977) is able to explain the ob- 
served extinction curves in the Milky Way as well as in the 
Local Group (e.g., Clayton et al. 2003). This model assumes 
uncoated refractory particles having a power law size distribu- 
tion n(a) - no a"^^, where a represents grain radius. As the 
present study is restricted to the optical regime, the model pro- 



posed by Mathis et al. (1977) serves as a good approximation 
for this purpose. 

Fig. |3] shows that R^ [= £(^717)] varies linearly with the 
inverse of the wavelength and is consistent with the result 
Qext °^ /I ' expected for small grains i.e., for x < \, where 
X - a 15 the grain radius and Q^xt is the extinction effi- 
ciency factor Thus, for a given Q^^, value, the mis-match be- 
tween the extinction curve for the Milky Way and those for the 
program galaxies is attributed to the difference in grain size be- 
tween the program galaxies and the Milky Way. Therefore, one 
can estimate the relative grain size in such cases by shifting the 
observed extinction curve along /I"' axis until it best matches 
the Galactic extinction curve (cf. Goudfrooij et al. 1994c), in 
the sense that the extinction curve lying below (above) the 
Galactic curve will correspond to smaller (larger) grain size 
relative to that of Galactic dust grains. The relative grain sizes 
thus derived for the sample galaxies are listed in Column 6 of 
TableS 

3.3. Dust mass estimation 

We have estimated dust content of the program galaxies using 
(i) optical extinction studied here and (ii) IRAS fluxes at 60 fim 
and 100 fim taken from Knapp et al. (1989), as described in the 
following subsections. 

3.3.1. Using optical extinction 

To estimate dust mass using total optical extinction, we fol- 
lowed the method described by Goudfrooij et al. (1994c),which 
is outlined here. For a given grain size distribution function 
n(fl) of the spherical grains of radius a and extinction efficiency 
Qextia, A), the extinction cross-section at wavelength A is given 
by 

Cext('i) — I Qext(o,A)7Ta n(d) da (2) 

where a_ and a+ are the lower and upper cutoffs of the grain 
size distribution, respectively. Assuming n{a) to be the same 
over the entire dusty region and using the definition of the ef- 
ficiency factor, Qext{a,A) - Cextia, A)/7Ta^ (ratio of extinction 
cross-section to the geometrical cross-section), the total extinc- 
tion due to dust at wavelength A is expressed as 

A, = 1.086 C,,,(^)xZrf (3) 

where Id is the dust column length along the line of sight. The 
column length density in units of g cm - for the dust is then 
expressed as 

'^d — I — 71 Pel n(a) da X 1^ (4) 

J a- ^ 

where gives the specific grain mass density which is taken 
to be ~ 3 g cm"^ for graphite and silicate grains (Draine & 
Lee 1984). This is then multiplied by the total area occupied 
by dust to obtain the dust mass, Md = 2d x Area, expressed in 
solar mass units. 
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Object 


Rb 


Rv 


Ri< 


Ri 


^?r> 

"Gal 

(6) 


(1) 


(2) 


(3) 


(4) 


(5) 


NGC 524 


4.46 ±0.39 


3.46 ±0.16 


2.21 ±0.47 


1.57 ±0.45 


1.06 ±0.02 


NGC 984 


4.19 ±0.36 


3.19 ±0.21 


2.43 ±0.19 


2.06 ±0.17 


1.08 ±0.08 


NGC 1172 


4.25 ±0.17 


3.25 ±0.19 


2.39 ±0.19 


1.95 ±0.14 


1.07 ±0.05 


NGC 1439 


3.45 ±0.21 


2.45 ±0.17 


2.01 ±0.15 


1.76±0.13 


0.93 ± 0.08 


NGC 2128 


4.37 ±0.28 


3.37 ± 0.27 


2.06 ± 0.29 


1.59 ±0.30 


1.04 ±0.06 


NGC 2534 


3.03 ±0.19 


2.03 ± 0.28 


1.61 ±0.30 


1.18±0.38 


0.80 ± 0.04 


NGC 2563 


— 


3.10 ±0.00 


2.29 ± 0.74 


1.19 ±0.51 


0.96 ± 0.03 


NGC 2672 


4.27 ±0.14 


3.27 ± 0.28 


2.48 ±0.14 


1.77 ±0.15 


1.06 ±0.02 


NGC 2693 


4.46 ± 0.23 


3.46 ± 0.26 


2.87 ± 0.45 


2.59 ±0.36 


1.30 ±0.04 


NGC 2907 


3.25 ±0.40 


2.25 ± 0.48 


1.19 ±0.29 


0.49 ± 0.23 


0.77 ± 0.03 


NGC 2911 


3.84 ±0.35 


2.84 ± 0.29 


1.98 ±0.31 


1.13 ±0.37 


0.92 ± 0.04 


NGC 3078 


4.10 ±0.17 


— 


2.39 ±0.19 


— 


1.02 ±0.03 


NGC 3489 


4.38 ±0.25 


3.38 ±0.21 


2.64 ±0.19 


1.70±0.15 


1.09 ± 0.02 


NGC 3497 




3 1 + no 




1.53 ± 0.65 


1.02 ± 0.02 


NGC 3585 


4.45 ± 0.44 


3.45 ± 0.33 


3.00 ± 0.24 




1.17 ±0.04 


NGC 3599 


3.49 ±0.15 


2.49 ± 0.09 


2.12 ±0.06 


1.57 ±0.14 


0.93 ± 0.06 


NGC 3665 


4.29 ± 0.20 


3.29 ±0.13 


2.33 ± 0.23 




1.03 ±0.02 


NGC 3923 


4.16 ±0.43 


3.16 ±0.39 


2.03 ±0.16 




0.98 ± 0.04 


NGC 4459 


3.94 ± 0.24 


2.94 ±0.18 


2.26 ±0.16 




0.90 ± 0.02 


NGC 5363 


3.78 ±0.27 


2.78 ±0.16 


2.18 ±0.36 


1.89 ±0.28 


0.92 ± 0.04 


NGC 5525 


4.15 ±0.29 


3.15 ±0.17 






0.99 ± 0.01 


NGC 5898 


4.15 ±0.30 


3.15 ±0.46 


2.46 ± 0.29 


1.84 ±0.41 


1.04 ±0.03 


NGC 5903 


4.07 ±0.30 


3.07 ± 0.22 


2.28 ± 0.39 


2.01 ±0.34 


1.05 ±0.08 


NGC 7432 


3.89 ±0.28 


2.89 ±0.21 


2.46 ±0.19 


1.94 ±0.22 


1.05 ±0.09 


NGC 7722 


3.89 ±0.25 


2.89 ± 0.20 


2.03 ±0.14 


1.78 ±0.10 


0.94 ± 0.02 






Rv 


Rr 


R; 




NGC 5485 


4.38 ± 0.24 3.65 ±0.21 


2.65 ±0.17 


1.86 ±0.12 


1.40 ±0.1 8 


0.92 ± 0.04 


The Galaxy 


4.10 


3.10 


2.27 


1.86 


1.00 



Notes on Table l4l Column (2) to (5) lists the R i ] values with A = B,V, R, &/. Column (6) lists the relative (mean) grain size 

in the program galaxies with respect to that responsible for the Galactic curve. NGC 5485 represents data in five bands {U,B,V,R & I). For 
comparison, Galactic R^ values are listed in the last row. 



The measured total extinction in V band (Ay) can be used 
to compute dust mass using the size distribution of Mathis et 
al. (1977) as 



n(a) — no a 



-3.5 



(a- < a < a+) 



where a_ = 0.005 /jm and a+ - 0.22 /j.m for the Milky Way 
with Ry - 3.1 (cf. Draine & Lee 1984). Since the observed 
extinction curves in the program galaxies refer to the dust 
grains at the upper end of the size distribution (Goudfrooij et 
al. 1994c), upper limits of the grain size in the target galaxies 
were scaled accordingly using the relation 



< a> 



flGal 



X 0.22 yUOT 



(5) 



where is the relative grain size for the program galaxies 
and are listed in Table |4] 

Further, to estimate the extinction efficiency of the dust grains 
we assume spherical grains composed of silicate and graphite 
with nearly equal abundance (see, Mathis et al. 1977). The val- 
ues of the extinction efficiency for silicate and graphite grains 
are taken as 



<2. 



ext, silicate — 



Qext,graphit 



0.^ a/asilicate fora < asilicate, 

0.8 for a > a^mcate 



2.0 a/a^i; 
2.0 



graphite 



for a < a 



graphites 



for a > Clgraphit 



with a silicate - 0.1 /vm, and Ographite - 0.05 yum. Using these 
parameters and the dust column density, the total dust mass 
(Md,opticai) contained in the program galaxies was estimated. In 
determining the mean visual extinction, we included all those 



regions with tv > 0.02. The computed dust masses for the sam- 
ple galaxies using total optical extinction (Md.opticai) are given 
in column 6 of Table |5] 

3.3.2. Using IRAS densities 

Using IRAS flux densities measured at 60 and 100 fim, we first 
calculate the dust grain temperature in the program galaxies 

using relation, Td = 49 d^f '"^ (Young et al. 1989). The dust 
content (MdjRAs) in the sample galaxies was then computed 
using the relation (Hildebrand 1983): 



4 n F 

Md = -apdD^ 

3 QMTd) 



(6) 



where a, pd and D are the grain radius, specific grain mass den- 
sity and distance of the galaxy, respectively. Fy, Qy & ByiJd) 
are the observed flux density, grain emissivity and the Planck 
function of temperature Td at frequency v, respectively. We 
adopted pd - 3 gcrrT^ and = 0.04 gcrrT^ for 0.1 jjm 
grains at 100 yum (Hildebrand 1983). The derived dust masses 
from IRAS flux densities are listed in column 5 of the Table 5. 
Due to the fact that IRAS was insensitive to the dust emitting 
at wavelengths longer than 100 fim (i.e., dust cooler than about 
20 K), these estimates of dust masses using IRAS flux densities 
represent the lower limits (Tsai & Mathews 1996). 

4. Discussion 

Extinction curves derived for the sample galaxies reveal that 
the Ry value lies in the range between 2.03 to 3.46 for el- 
lipticals, and 2.25 to 3.46 for lenticulars, with an average 
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Fig. 3. Extinction curves for the program galaxies (filled squares, solid lines) along with the canonical curve for the Galaxy (open 
circles, dashed lines) for comparison. The offset (shift of the curves along the X or Y axis in order to separate them) for the 
X-axis is 4, that for the Y-axis is 3. 



of 3.01, compared to 3.1 in the case of the Milky Way. 
For some of the sample galaxies the Ry value is close to 
the classical Galactic value, while for some of them the Ry 
value is significantly less than 3.1. Interestingly, all galaxies 
(e.g., NGC 2534, 2907, 291 1, 5363, 5485 and 7722) show- 
ing smaller Ry values exhibit well-settled dust lane or lanes 
with an average Ry equal to 2.80, implying that the "larger" 
grains responsible for the optical extinction are roughly 
25% "smaller" in size compared to the grains responsible 
for the extinction curve in the Milky Way. For a sample of 
10 elliptical galaxies, Goudfrooij et al. (1994c) arrived at 
Ry in the range of 2.1 to 3.3, with characteristic dust grains 
up to 30% smaller in size compared to the canonical grains 
in the Milky Way. Thus, our results are in good agree- 
ment with those of Goudfrooij et al. (1994c). For a few of 
the remaining galaxies (e.g., NGC 524, 2672, 2693, 3489, 
3585) extinction curves are found to lie slightly above the 
Galactic extinction curve, implying marginally larger val- 
ues of Ry, and these are found to have ring or arc shaped 
dust morphologies. Thus, our results seems to indicate that 
the galaxies having Ry values smaller than the canonical 
value (and hence smaller grains) exhibit smooth, regularly 
distributed dust lanes, whereas the galaxies with larger Ry 



values exhibit irregular dust morphologies (Goudfrooij et 
al. 1994c). 

Dust masses derived from optical extinction are found to lie 
in the range ~ 10"* to 10^ M0 and are in good agree- 
ment with the earlier estimates for the early-type galaxies 
(Brosch et al. 1990; Goudfrooij et al. 1994c; Sahu et al. 
1998; Ferrari et al. 1999; Dewangan et al. 1999; Tran et 
al. 2001; Patil et al. 2002). The dust masses derived using 
the optical methods depend on the assumption of a fore- 
ground screen in front of the background light. Therefore, 
this method is not sensitive to any dust embedded or inter- 
mixed with the stars within the galaxy and hence the opti- 
cal method always provides a lower limit of the true dust 
content of the host galaxy. Comparison of the dust masses 
derived from IRAS flux densities and optical extinctions 
reveals that MdjRAS is larger than Md,opticd by a factor of 
7.8 for the sample of galaxies studied here. Goudfrooij & 
de Jong (1995) for a sample of elliptical galaxies alone re- 
ported this ratio to be equal to 8.4. In the case of SO galax- 
ies, as pointed out by Sahu et al. (1998) and Dewangan et al. 
(1999), this mass discrepancy is less significant compared 
to that for ellipticals. Using ISO data Temi et al. (2004) 
estimated a dust content of NGC 1172 and NGC 5363 of 



Table 5. Dust properties 
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Object 


dust moiphology 


IRAS flux (mJy) 


Log(M J.IRAS ) 


Log(Mj,„pii„,i) 


Td 


M 








60 ± err 


100 ± eiT 


(Mq) 


(Mq) 


(K) 


(Moyr-') 


(10" yr) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


NGC 524 


dust rings 


780±34 


1820±114 


5.95±0.11 


4.91 ±0.1 4 


35±4 


1.06 


7.8 


NGC 984 


maj. dust lane 


140±37 


I20±106 


4.85±0.88 


4.50±0.07 


< 30 


<0.89 


5.5 


NGC 1172 


nucl. dust 


0±41 


0±72 


<4.35 


4.16±0.18 


- 


< 0.12 


13.5 


NGC 1439 


min. dust lane 


0± 32 


300±35 


< 5.80 


5.43±0.19 


<20 


<0.04 


23.0 


NGC 2128 


min.dust lane 


870±37 


2160±166 


6.22±0.11 


4.42±0.12 


34±3 


0.10 


- 


NGC 2534 


min. dust lane 


380±50 


800±184 


5.71±0.28 


4.10±0.18 


36±4 


0.15 


- 


NGC 2563 


min. dust lane 


0±27 


± 187 


<4.70 


3.75±0.15 


- 


<0.51 


3.3 


NGC 2672 


nucl. dust 


0±45 


440±45 


< 6.30 


5.30±0.30 


<20 


< 0.69 


- 


NGC 2693 


min dust lane 


210±37 


790±114 


6.09±0.20 


4.84±0.15 


29±4 


0.02 


5.8 


NGC 2907 


multi. lanes 


310±26 


1090±194 


5.70±0.21 


5.05±0.22 


30±2 


0.47 


- 


NGC 2911 


maj. dust lane 


290±28 


560±79 


5.59±0.19 


5.54±0.14 


38±5 


0.10 


8.8 


NGC 3078 


nucl. dust 








3.91±0.08 






11.0 


NGC 3489 


dust clouds 








3.94±0.26 






18.5 


NGC 3497 


multi. lanes 


278 


< 1347 


< 6.49 


5.08±0.14 


< 26 






NGC 3585 


nucl. dust 


160±42 


0±70 


<4.24 


3.73±0.25 


< 30 


0.37 


14.7 


NGC 3599 


min. dust lane 








3.73±0.23 






24.0 


NGC 3665 


maj. dust lane 


1960±40 


6690±163 


6.60±0.07 


5.02±0.28 


30±3 


0.19 


23.1 


NGC 3923 


nucl. dust 


0±35 


0± 120 


<4.70 


4.46±0.08 




0.61 


10.6 


NGC 4459 


dust ring 


1920±67 


4280±119 


5.71±0.12 


3.66±0.14 


36±8 


0.48 


23.2 


NGC 5363 


multi. lanes 


1700±46 


4450±45 


5.72±0.18 


5.05±0.14 


33±13 


0.01 


23.1 


NGC 5485 


min. dust lane 


150±34 


850±88 


6.01±0.17 


4.36±0.11 


25±4 


0.16 


13.0 


NGC 5525 


maj. lane + ring 








5.99±0.10 








NGC 5898 


min dust lane 


130±37 


200±64 


4.85±0.41 


4.48±0.07 


41±9 


0.25 


5.4 


NGC 5903 


nucl. dust 


± 23 


0± 112 


< 6.31 


4.05±0.08 




0.43 


5.2 


NGC 7432 


min. dust lane 








3.98±0.09 






5.8 


NGC 7722 


multi. lanes 


820±34 


2840±159 


6.82±0.08 


6.19±0.16 


30±2 


1.72 





Notes to Table [s] Column(2) lists morphology of dust: dust ring(s) - face on dust ring(s); maj. dust lane - dust lane along the optical major axis; min. 
dust lane - dust lane aligned along minor axis; multi. lanes - multiple dust lanes; dust clouds - dust distribution in the form of clouds; nucl. dust - nuclear 
dust. Column (3) & (4) list IRAS flux densities at 60/jm and 100/im (taken from Knapp et al. 1989). Column (5) and (6) list dust masses derived from 
IRAS flux densities (Mj iras) and optical extinction values (Md.opiioai), respectively. Column (7) gives dust temperature (Tj), while column (8) lists the 
present day mass-loss rate from the evolved, red-giant stars in the program galaxies. Column (9) lists life time of dust grains derived using Eg. ([8]. 



9.9 X 10^ Mq & 2.00 X 10^ Mq, respectively, these estimates 
are roughly an order of magnitude larger than our estimates 
from IRAS fluxes. This in turn implies that although ISO 
data acts as a more reliable tracer of the true dust content in 
the extragalactic environment, it renders the dust mass dis- 
crepancy even more significant than using IRAS flux den- 
sities. 

There are several processes that may alter the effective dust 
grain size in the external galaxies, like destruction of grains 
due to sputtering in supernova blast waves, grain-grain col- 
lisions, sputtering by thermal ions (warm and hot) etc., 
as has been discussed by Goudfrooij (1999). But, which 
one of these would work as the dominant process for de- 
struction of grains leading to smaller grain size depends on 
the conditions of specific galaxies. In particular, for galax- 
ies immersed in hot X-ray emitting gas, sputtering by hot 
ions is expected to be the dominant process for destruction 
and alteration of the grain size. In order to investigate this 
point further we collected X-ray flux for the sample galax- 
ies from the literature, and examined the relationship be- 
tween the dust mass and X-ray luminosity, both normalized 
to the blue luminosity. The result shown in Fig.SJa) reveals 
a marginal anti-correlation suggesting that high (Lx/Lb) 
galaxies harbour less dust, in agreement with the predic- 
tions of Goudfrooij et al. (1994c). 

The issue of the origin of dust in early-type galaxies is rather 
controversial. Suggested sources by which galaxies may 
acquire dust include (i) internal: in the form of mass-loss 
from evolved stars distributed throughout the galaxy more 
or less uniformly, and (ii) external: in which a gas-rich com- 



panion donates gas and dust either via direct merger, accre- 
tion or tidal capture. Dwek & Scalo (1980) have demon- 
strated that evolved red giant stars are the dominant source 
of dust injection into the ISM of the early-type galaxies. To 
check the significance of internal sources, using the present 
day mass-loss rate we estimated the total amount of dust 
accumulated over the lifetime of a galaxy under the as- 
sumption that the properties of the evolved, red-giant stars 
in the Galactic bulge closely resemble those in elliptical 
galaxies. For most of the K & early M giant stars with ef- 
fective surface temperatures about 4500 K and with no cir- 
cumstellar material, one finds S n ^im IS 2.2 - 0.08 (Knapp 
et al. 1992), while for the stars with circumstellar matter (as 
in elliptical galaxies) this ratio is 0.13, about 38% larger 
than expected from stars in the Galactic bulge (Knapp et 
al. 1992). This excess flux measured at 12 jjm in elliptical 
galaxies is attributed to the emission from the circumstellar 
shells produced by mass-loss from the evolved stars, and 
thus provides a reasonable estimate for the total mass-loss. 
Following Knapp et al. (1992), we adopt the following re- 
lation for the mass-loss rate: 



M = 6xl0" 



— y 

Mpc j 



5,2 -0.042 5 



100 



mJy 



MQyr- 



(7) 



The term 0.042 Sioo gives the correction to the excess flux 
at 12 /vm emission due to contamination from the interstel- 
lar dust heated to an equilibrium temperature of about 20 K 
by the ambient ultraviolet interstellar radiation field. The 
inferred values of the mass-loss rates in the sample galax- 
ies are listed in column 8 of Table |5] 
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Fig. 4. The ratio of dust mass to blue luminosity (— as a function of (a) ratio of the X-ray-to-blue luminosity, and (b) blue 
luminosity (filled squares are for E, while open triangles for SO). 



Dust injected into the ISM by the evolved stars through 
mass-loss is also destroyed simultaneously by a variety of 
processes. In the absence of hot. X-ray emitting gas, the 
life time of the 0. 1 jum refractory grains against sputtering 
through low velocity shock waves or grain-grain collision 
is ~ 10'^ yr, while that in the presence of hot coronal gas 
(Te ~ 10^ K) against collision with the hot protons and a- 
particles is of the order of ~ lO*" - 10^ yr (cf. Draine & 
Salpeter 1979, Goudfrooij & de Jong 1995). A majority of 
the galaxies from the present sample are detected as X-ray 
emitters, and therefore, we estimate the lifetime of the dust 
grains (column 9 of Table |5]l of radius a against such colH- 
sion using the relation 



2x lO'* 



^-3 



yr 



(8) 



(Draine & Salpeter 1979), where hh is the proton density 
in the plasma and a is the grain size. Here rin - 0.83 rie, rie 
being the electron density estimated using 



n^ir) = n^(0) 
and 

nM = 0.061 



1+1^ 



-3/4 



10'*'ergs s 



1/2 



1 kpc 



-3/2 



cm 



-3 



where rieiO) is the central number density over the X-ray 

core radius ax(^- 1-0 ("[oTt^) kpc) (see, e.g., Canizares 
et al. 1987, for details). X-ray luminosities reported by 
Bettoni et al. (2003) were used to derive these parameters. 
For the grains of radius 0. 1 fim this time scale is found to 
lie in the range between 2.4 x lO^yr and 3.3 x 10^ yr which 
correspond to NGC 3599 and NGC 2563 having the lowest 
and highest X-ray luminosities, respectively, in the present 
sample. 



For illustration, we examine the accumulation of dust for 
the case of NGC 5363 (log Lx= 40.38 ergs s^', nn = 2.2 x 
lO-\m-\Td = 2.3 X 10^ yr and M = lO^^Moyr"' and 
gas-to-dust ratio of about 100), taking into account the two 
competing processes, namely mass-loss in the form of dust 
from the evolved stars and destruction of dust. The rate of 
accumulation of dust is 



dt 



where 



dt 



Md(t) T-. 



is the rate of mass loss in the form of dust from 
red giant stars and r^' is the destruction rate. Integrating 
this dust accumulation rate over time, the total build up of 
dust mass in NGC 5363 over its life time is shown in Fig.|5] 
and is found to be equal to 8.5 xlO-'M©, which is about 
a factor of 100 times smaller than the observed dust mass. 



10000 



8000 



6000 



4000 



EOOO 




10= 2x108 3x10" 

Time (yrs) 

Fig. 5. Total build up of dust mass within NGC 5363 over its 
life time 
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For other galaxies, where X-ray flux was available in the lit- 
erature, comparison of the build up mass and the observed 
mass shows a similar discrepancy except for NGC 3585 and 
NGC 3923, for which the two estimates are roughly equal. 
This in turn implies that internal sources are not sufficient 
to account for the observed dust and therefore an external 
source plays a dominant role in the supply of dust into the 
ISM of early-type galaxies. 
There are several other indicators that support the external 
origin of dust in the early-type galaxies. If it is of internal 
origin, it should exhibit the same global rotational dynam- 
ics as that of the stellar system (Kley & Mathews 1995). 
However, kinematical studies of dust lane galaxies have 
shown that the angular momentum vectors of the inter- 
stellar gas are often orthogonal to those of the stars, im- 
plying that the ISM must have a different origin (Caon et 
al. 2000). Another good argument in favour of an exter- 
nal, merger-related origin of dust in early-type galaxies is 
the lack of significant correlation between the Lfm or dust 
mass (Md iRAs) and optical luminosity {Lb), Fig.HJb). This 
suggests that at least some of our sample galaxies have ac- 
quired dust externally in a merger-like event. Similar con- 
clusions were also arrived at by Forbes (1991), Goudfrooij 
& de Jong (1995), Temi et al. (2004). 




6.5 7 

Fig. 6. Environmental dependence of dust mass; upper and 
lower panels represent dependence on group members (N) and 
group harmonic radius (Rh), respectively 

To explore the possibility that dust could originate through a 
direct merger or accretion or tidal capture of gas from their 
encounter with neighbouring galaxies, we searched the en- 
vironment of the program galaxies and found that the ma- 
jority of the galaxies are either a part of a group / cluster, or 
form interacting pairs with their neighbours. Relevant de- 
tails on the environment of the galaxies taken from Huchra 
& Geller (1982), Gefler & Huchra (1983) and Gai'cia 
(1993) are given in the column 8 of Table[T] Following Tran 
et al. (2001), we also examined the relationship between the 
dust content and environment of the program galaxies using 
parameters like group membership (A^) and group harmonic 
radius (Rh)', Rh provides a rough measure of the size of the 



group (see, Tran et al. 2001 for details). The group mem- 
bership and group harmonic radius is plotted as a function 
of dust mass in Fig. |6] The scatter plots between vs. 

& Rh indicate the absence of any significant correlation 
between dust mass and the two quantities. As group mem- 
bership in these catalogs mostly relies on the radial veloc- 
ity measurement of high-luminosity galaxies, small (dwarf) 
galaxies are often missing from these catalogs, and hence 
log N may not be representative of the real situation. If the 
donors of the dust are supposed to be galaxies similar to 
those present in the group catalogs mentioned above, this 
would make the possibility of dust originating from neigh- 
bouring companions through merger or tidal capture more 
complicated. However, one should interpret these results 
with caution, as a given dusty galaxy might have acquired 
dust sufficiently long ago that no signature of merger or in- 
teraction is visible at the present epoch (Tran et al. 2001). 

5. Conclusions 

We have reported properties of dust in a sample of 26 early- 
type galaxies based on their multicolour imaging observations. 
Our results are summarized here. 

Extinction curves derived for galaxies studied here run par- 
allel to the canonical curve of the Milky Way, implying that 
properties of dust in the extragalactic environment are similar 
to those of the canonical grains in the Milky Way. The Ry value, 
which characterizes the extinction curve in the optical region, 
is found to vary in the range 2.03 - 3.46 with an average of 
3.02, and is not very much different from the canonical value 
of 3.1 for our Galaxy. However, several galaxies studied here 
exhibit smaller Ry values suggesting that the "large" grains re- 
sponsible for the optical extinction are significantly smaller in 
size than the canonical grains in the Milky Way. Our results in- 
dicate that Ry values largely depend on the morphology of dust 
in the host galaxies, in the sense that galaxies with well-defined 
dust lanes exhibit much smaller grains compared to those with 
irregular dust morphologies. 

The dust content derived by the optical extinction method 
for the program galaxies is found to lie in the range 10"* to 
10^ M0 and is, in general, less than the dust mass derived from 
IRAS flux densities. This discrepancy is more significant in the 
case of ellipticals than lenticulars, in good agreement with the 
results obtained by other researchers for early-type galaxies. 

Galaxies with well settled dust lanes contain high dust 
masses and are weak X-ray sources, in agreement with the 
prediction of Goudfrooij (1994c). Further, galaxies with well- 
settled dust lanes have smaller values of Ry, and are found to 
lie in denser environments. 

For origin of dust in early-type galaxies, internal mass-loss 
from evolved stars alone cannot account for the observed dust 
in these galaxies and therefore points towards an external origin 
for the dust in these galaxies. 
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